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This work presents the first continuum shell-model study of weakly bound neutron-rich nuclei involv- 
ing multiconfiguration mixing. For the single-particle basis, the complex-energy Berggren ensemble 
representing the bound single-particle states, narrow resonances, and the non-resonant continuum 
background is taken. Our shell-model Hamiltonian consists of a one-body finite potential and a 
zero-range residual two-body interaction. The systems with two valence neutrons are considered. 
The Gamow shell model, which is a straightforward extension of the traditional shell model, is shown 
to be an excellent tool for the microscopic description of weakly bound systems. It is demonstrated 
that the residual interaction coupling to the particle continuum is important; in some cases, it can 
give rise to the binding of a nucleus. 
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The microscopic structure of exotic nuclei near the 
particle drip lines is a topic of great current interest in 
low-energy nuclear physics. Apart from theoretical and 
experimental nuclear structure interest, calculations for 
nuclei far from stability have astrophysical implications, 
especially in the context of stellar nucleosynthesis. What 
makes this subject particularly challenging is the weak 
binding; hence the closeness of the particle continuum. 

There are many factors which make the coupling to 
the particle continuum important. Firstly, even for a 
bound nucleus, there appears a virtual scattering into 
the phase space of unbound states. Although this pro- 
cess involves intermediate scattering states, the corre- 
lated bound states must be particle stable, i.e., they have 
zero width. Secondly, the properties of unbound states, 
i.e., above the particle (or cluster) threshold dirccly re- 
flect the continuum structure. In addition, continuum 
coupling directly affects the effective nucleon-nucleon in- 
teraction. 

The treatment of continuum states is an old problem 
which, in the context of excited states near or above the 
decay threshold, has been a playground of the contin- 
uum shell model (CSM) §. In the CSM, including the 
recently developed Shell Model Embedded in the Contin- 
uum (SMEC) the scattering states and bound states 
are treated on an equal footing. So far, most applications 
of the CSM, including SMEC, have been used to describe 
situations in which there is coupling to one-nucleon decay 
channels. However, by allowing only one particle to be 
present in the continuum, it is impossible to apply the 
CSM to 'Borromean systems' for which A- and (^-2)- 
nucleon systems are particle-stable but the intermediate 
(A-l)-system is not. Several approaches, including the 
hyperspherical harmonic method or the coupled-channel 



approach, have been developed to study structure and 
reaction aspects of three-body weakly bound nuclei [Q. 
However, most of these models utilize the particle-core 
coupling which does not allow for the exact treatment of 
core excitations and the antisymmctrization between the 
core nucleons and the valence particles. 

The reason for limiting oneself to only one particle in 
the continuum is two-fold. First, the number of scat- 
tering states needed to properly describe the underlying 
dynamics can easily go beyond the limit of what present 
computers can handle. Second, treating the continuum- 
continuum coupling, which is always present when two 
or more particles are scattered to unbound levels, is dif- 
ficult. There have been only a few attempts to treat the 
multi-particle case and, unfortunately, the proposed 
numerical schemes, due to their complexity, have never 
been adopted in microscopic calculations involving mul- 
ticonfiguration mixing. Consequently, an entirely differ- 
ent approach is called for. In this work, we formulate 
and test the multiconfigurational shell model in the com- 
plete Berggren basis. The resulting Gamow Shell Model 
(GSM) is then applied to systems with two valence neu- 
trons. 

The Gamow states (sometimes called Siegert or reso- 
nant states) are generalized eigenstates of the time- 
independent Schrodinger equation with complex energy 
eigenvalues E = Eq — iT/2, where F stands for the de- 
cay width (which is zero for bound states). These states 
correspond to the poles of the S'-matrix in the complex 
energy plane lying on or below the positive real axis; they 
are regular at the origin and satisfy a purely outgoing 
asymptotics. In the following, we consider the Gamow 
states of a one-body spherical finite potential. The single- 
particle (s.p.) basis of Gamow states must be completed 
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by means of a set of non-resonant continuum states. This 
completeness relation, introduced by Berggren 0, reads: 



|</.,(fc))(</>,(r)|dfc = i, (1) 
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where (^nj are the Gamow states carrying the s.p. angu- 
lar momentum n stands for all the remaining quantum 
numbers labeling Gamow states, 4>j{k) are the modified 
scattering Gamow states, and the contour i+ in the com- 
plex fc-plane has to be chosen in such a way that all the 
poles in the discrete sum in Eq. are contained in the 
domain between and the real energy axis. If Unj{r) 



stands for the radial part of (pnj, then Unj{r) 



dry 



and (j)nj—(f>nj{u — > m). If the contour i+ is chosen rea- 
sonably close to the real energy axis, the first term in 
(|l|) represents the contribution from bound states and 
narrow resonances while the integral part accounts for 
the non-resonant continuum. A number of completeness 
relations similar to (|l|) were studied by Lind Q]. 

There have been several applications of resonant states 
to problems involving continuum j^, but in most cases 
the so-called pole expansion, neglecting the contour inte- 
gral in Eq. (|l ) , was used . The importance of the con- 
tour contribution was investigated in Rcfs. in the 
context of the continuum RPA with separable particle- 
hole interactions where it was concluded that the non- 
resonant part must be accounted for if one aims at a 
quantitative description. This can be achieved by dis- 
cretizing the integral in Eq. ([|) |l2| : 

N 

|0,(fc))(0,(r)|dfc = ^l<^,(fc,))(<^,(fc*)|A,, (2) 

where depends on the quadrature used (in our case 
we use the four-point interpolation). 

In our study, Gamow states are determined using the 
generalized shooting method for bound states which re- 
quires an exterior complex scaling The numerical 
algorithm for finding Gamow states for any finite-depth 
potential U{r) has been tested on the example of the 
Poschl-Teller-Ginocchio (PTG) potential for which 
the resonance energies and wave functions are known an- 
alytically. Energies of all PTG resonances with a width 
of up to 90 MeV are reproduced with a precision of at 
least 10^^ MeV. The antisymmetric two-particle wave 
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)j are obtained in the usual way by 
coupling the s.p. wave functions of the considered bound, 
resonance, and scattering Gamow states labeled by sub- 
scripts ii, i2 to the total angular momentum J. The com- 
pleteness relation for two-particle states. 
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can be used to calculate the two-body matrix elements. 
The radial integrals entering the Hamiltonian matrix el- 
ements were regularized separately by an appropriate 



choice of the angle of the external complex scaling. The 
resulting (complex symmetrix) Hamiltonian matrix can 
be diagonalized using standard methods. 

In most applications, one is interested in bound or reso- 
nance A^-body states but not in non-resonant continuum. 
Bound states can be clearly identified, because the imagi- 
nary part of their energy must be zero. No equally simple 
criterion exists for resonance or scattering states. On the 
other hand, the coupling between scattering states and 
resonant states is usually weak, so one can determine the 
resonances using the following two-step procedure. In the 
first step, the shell-model Hamiltonian is diagonalized in 
both (i) the full space including the contour, and (ii) the 
subspace of Gamow states (pole expansion). In the sec- 
ond step, one identifies the eigenstates of (i) which have 
the largest overlap with those of the second diagonaliza- 
tion. For the case of two valence particles discussed in 
this work, one can include in the basis up to 50 states 
in the non-resonant scattering continuum. For greater 
dimensions, e.g., for a larger number of valence particles, 
this method becomes impractical and the perturbative 
correction methods must be used [ pT[ . 

In the following exploratory GSM calculations, we shall 
consider two cases: (i) ^^O with the inert -'^^O core and 
two active neutrons in the sd shell, and (ii) ^He with the 
inert "^He core and two active neutrons in the p shell. Our 
aim is not to give the precise description of ^^O and ^He 
(for this, one would need a realistic Hamiltonian and a 
larger configuration space), but rather to illustrate the 
method, its basic ingredients, and underlying features. 

The "1^0" case 

The s.p. basis was generated by a Woods-Saxon (WS) 
potential with the radius i?o=3.05fm, the surface diffuse- 
ness (i=0.65fm, the potential depth [/o=-55.8MeV, and 
the strength of the spin-orbit term [7^0=6.06 MeV. With 
this choice of parameters, the single particle 0^5/2 and 
lsi/2 states are bound, with s.p. energies -4.14MeV, and 
-3.27 MeV, respectively, and 0c?3/2 is a resonance with the 
s.p. energy 0.9-10.97 MeV. Energies of these s.p. states 
are close to the s.p. states of ^''O. 

The completeness relation requires taking the Si/2, 
c?5/2, and ^3/2 non-resonant continuums. For the lsi/2 



and bound states, their non-resonant continuums 

can be chosen along the real momentum axis. Since, to 
the first order, the inclusion of these continuums should 
only result in the renormalization of the effective inter- 
action, they are ignored for the purpose of the present 
exercise whose main focus is the neutron emission. On 
the contrary, 0^3/2 is a resonance state, so the associated 
contour has to be complex to produce the correct energy 
width. The contour representing the (i3/2-continuum 
was chosen to consist of three straight segments con- 
necting the points ki~Q-\Q, fc2=0.3-0.2i, A:3=0.5-i0, and 
/c4=2.0-i0 (all in fm"^). The strength of the (5-force was 
taken to be Vb=-350MeVfm3. 
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The completeness of the Gamow basis depends on the 
number of discretized scattering basis states considered. 
Table | illustrates this dependence. The real part of en- 
ergy represents the binding energy of a state with respect 
to the ^^O core, i.e., the two-neutron separation energy. 
For the resonance states, the real and imaginary parts 
of energy do not change much by increasing the number 
of scattering states. On the other hand, bound states 
acquire a very small negative width which does not ex- 
ceed several keV. This spurious negative width depends 
strongly on the basis size, and the convergence to zero is 
both slow and non-monotonic. The presence of a small 
and negative width is a feature of particle-bound states 
obtained in the GSM. The results displayed in Table || 
show that only about 10-20 vectors in the scattering con- 
tinuum are sufficient to keep the error of calculated en- 
ergies and widths at the acceptable level. It is also clear 
that the "no-pole" approximation (inclusion of no scat- 
tering states) gives a rather poor description of bound 
and near-threshold states while it works fairly well for 
high-lying states carrying a sizeable width. In this re- 
spect, this result is consistent with the conclusions of 
Refs. [0Jll). 

In the considered example, the calculated one-neutron 
threshold is -4.142 MeV (0^5/2 energy) while the two- 
neutron threshold is at zero (the binding energy of the 
core). Consequently, few states shown in Table | are 
two-neutron bound but one-neutron unstable, e.g. the 
calculated width, 57keV, of the J'^=2^ state at i?o=- 
2.61 MeV, characterizes single neutron emission from this 
state. The higher-lying states shown in Table | are unsta- 
ble with respect to both one- and two-neutron emission. 

TABLE I. Dependence of energies (left number, in MeV) 
and neutron widths (right number, in keV) of calculated 
states in ^**0 on the number of discretized scattering basis 
states along the contour L+. 



culated for the bound {Of and 2^ ) and resonance {2^) 
states in ^^O. All eigenstates are normalized according to 
Berggren (|J^: J2n = 1- One should notice that - con- 
trary to the conventional SM case ~ no modulus square 
appears in the normalization. This implies that the prob- 
abilistic interpretation of c„ must be generalized , 
i.e, when computing expectation values the real part of 
should be associated with the mean value while the imag- 
inary part represents the uncertainty due to the decaying 
nature of the state. As seen in Table |l|, the contribution 
from the non-resonant continuum plays a different role 
compared to that from the resonant states. Firstly, it 
is generally smaller than the leading components involv- 
ing resonant orbits, though in the example shown in Ta- 
ble |l| the contribution of the 0dy2 resonance is similar in 
magnitude. Secondly, according to our calculations, the 
inclusion of the contour primarily affects the imaginary 
part. Finally, the contribution from two particles in the 
non-resonant contimmm, L:^_ \ even though smaller than 
the one-particle contribution, L^^\ is not negligible. 

TABLE II. Squared amplitudes of different configurations 
in 0^, 2^, and 2f states of ^^0. The sum of squared ampli- 
tudes of all Slater determinants including one and two parti- 
cles in the non-resonant continuum are denoted by L^^^ and 
(2) 

L\ , respectively. 50 discretized scattering states were used. 
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■'5/2 
0^3/2 
lSl/20d5/2 
lSl/20d3/2 
04/20^3/2 



L 



(2) 



or 



2^ 



0.05-i9.1E-6 
0.91-i6.1E-6 0.86-Hil.2E-5 
0.02-i5.3E-3 1.9E-3-i4.4E-4 
0.13-il.2E-5 
4.6E-3-i5E-4 
7.7E-3-i8.4E-4 
1.3E-2-fi3.8E-3 2.5E-3-hil.7E-3 
3.4E-3+il.5E-3 1.6E-4-hi9.9E-5 



6.9E-3-i5.3E04 
1.6E-3-i5.2E-4 
4.5E-3-i3.6E-4 
0.03-i4.9E-3 
0.96-hi4.5E-3 
-1.2E-3-Hil.7E-3 
5.5E-5-hi4.3E-5 



rj TV 


states 


10 states 


30 states 


50 states 






Bound States 




ot 


-11.73, -131 


-12.11, -2.91 


-12.12, 0.27 


-12.12, 0.21 




-9.20, -26.37 


-9.24, -0.51 


-9.24, -0.031 


-9.24, -0.032 




-8.64, -13.51 


-8.64, -0.25 


-8.64, -0.004 


-8.64, -4E-4 




-7.66, -1.08 


-7.66, -0.324 


-7.66, -0.264 


-7.66, -0.260 


2j 


-7.85, -4.64 


-7.86, -0.167 


-7.86, -0.066 


-7.86, -0.049 






Resonances 




4J 


-4.00, -89.8 


-4.03, -8.15 


-4.04, 0.54 


-4.04, 0.54 


2+ 
^3 


-3.48, 29.08 


-3.48, 43.08 


-3.49, 44.65 


-3.49, 44.65 


24+ 


-2.60, 46.45 


-2.61, 55.82 


-2.61, 57.02 


-2.61, 57.04 




1.05, -172 


0.94, -18.64 


0.94, 0.82 


0.94, 0.76 


n 


1.63, 112.4 


1.64, 111.9 


1.63, 117.5 


1.63, 118.5 



In order to illustrate the configuration mixing in- 
duced by the two-body interaction in the GSM, Table || 
shows the complex squared shell-model amplitudes cal- 



The "^He" case 

A description of the Borromean nucleus ^He is a challenge 
for the GSM. "^He is a well-bound system with the one- 
neutron emission threshold at 20.58 MeV. On the con- 
trary, the nucleus ^He, with one neutron in the p shell, is 
unstable with respect to the neutron emission. Indeed, 
the J"" = ground state of ^He lies 890 keV above 

the neutron emission threshold and its neutron width is 
large, F^GOOkeV. The first excited state, 1/2^, is a very 
broad resonance (F=:4 MeV) that lies 4.89 MeV above the 
threshold. ^He, on the contrary, is bound with the two- 
neutron emission threshold at 0.98 MeV and one-neutron 
emission threshold at 1.87 MeV. The first excited state 
2'1 at 1.8 MeV in ^He is neutron unstable with a width 
F=113keV. In our GSM calculations, the states in ^He 
are viewed as one-neutron resonances outside of the ''He 
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core. A good fit to 3/2;^ and 1/2^ states in ^He is ob- 
tained by taking the WS potential with Rq = 2.0 fm, 
d=0.65fm, Uo = -47.0 MeV, and Uso = 7.5 MeV. With 
this potential, one finds the single-neutron resonances 
P3/2 and pi/2 at i;=0.745-i0.32 MeV and i;=2. 13-12.94 
MeV, respectively. The s.p. basis has been restricted to 
the Op3/2 resonance state and the P3/2 non-resonant con- 
tinuum. [The Opi/2 resonance is very broad and cannot 
be included in a meaningful way in the discrete sum in 
Eqs. (0,0). Consequently, following the reasoning applied 
to the^ case, the pi/2 contour along the real fc-axis has 
been ignored.] The L+-contour for the non-resonant P3/2 
continuum is chosen to enclose the 0^3/2 resonance: fci = 
O-iO, k2 = 0.2-0.2i, fc3=0.5-i0, and /c4=2.0-i0 fm'^ The 
strength of the (5-force was taken to be Vb=650MeVfm^. 
The number of points used to discrctize the scattering 
continuum is 50, though even with 15 points the results 
are reasonably stable. With this precision, we reproduce 
the most important feature of ®He: the ground state is 
particle-bound, despite the fact that all the basis states lie 
in the continuum. Table III shows the structure of wave 
functions of 0^ and 2f states in ^He. The important 
contribution from the non-resonant continuum is seen, 
even for the 0^ ground state which is particle stable. In 
spite of a very crude Hamiltonian, the neglect of the exact 
three-body asymptotics, etc., the calculated ground state 
energy £'=-0.951-10.01 MeV reproduces surprisingly well 
the experimental ground state energy with respect to the 
two-neutron emission threshold. The excited state 2f is 
predicted to lie at 2.25 MeV, slightly above the exper- 
imental value, and its width r=700keV which depends 
sensitively on the position of the state with respect to 
the emission threshold is somewhat too high as well. 

TABLE III. Same as in Table 0, except for 0+ and 2+ 
states in ''He. 







2^ 


OP3/2 


0.95-10.79 


l.Oll-fiO.0044 




O.ll+iO.76 


-0.011-10.0049 




-0.06+i0.03 


-1.3E-4-Hi4.8E-4 



In conclusion, the complex-energy Berggren ensemble 
is applied for the first time in shell-model calculations 
for two-neutron states near the particle-emission thresh- 
old. The results are very encouraging. It is seen that the 
contribution from the non-resonant continuum is impor- 
tant, especially for bound and near-threshold states. The 
particle-bound states calculated in the GSM are charac- 
terized by small and negative widths which show non- 
monotonic behavior as a function of the basis size. Ac- 
cording to our experience, only about 10-20 vectors in 
the scattering continuum are sufficient to keep the er- 
ror of calculated energies and widths at an acceptable 



level. With a simple interaction, such as the 5-force, 
we calculated the low- lying states of and ^He and 
discussed their properties with respect to neutron emis- 
sion. Last, but not least, pairing correlations due to 
the continuum-continuum scattering have been shown to 
bind the ground state of ^He with a completely unbound 
basis provided by the s.p. resonances of ^He. Further 
applications of the GSM are in progress [H . 
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